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ABSTRACT 


The  paper  presents  the  results  of  new,  precise  measure¬ 
ments  of  the  viscosity  of  nitrogen,  argon,  and  helium  at  25°C. 
The  measurements  were  performed  over  a  nominal  range  of  pres¬ 
sures  1-100  atm  and  at  very  closely  spaced  density  intervals. 
The  data  were  subjected  to  a  stringent  statistical  analysis 
in  order  to  determine  whether  the  density  expansion  consists 
of  a  pure  polynomial  or  whether  a  term  of  the  form  p2  Jtnp 
must  be  included  in  it.  The  existence  of  such  a  term  was  dis¬ 
covered  theoretically  by  several  investigators.  The  analysis 
indicates  that  if  such  a  term  exists,  its  factor  must  be  very 
small . 

This  result  is  interpreted  as  indicating  that  correlations 
which  extend  over  distances  of  the  order  of  a  mean  free  path 
are  negligible  when  compared  with  correlations  which  extend 
over  distances  of  the  order  of  the  range  of  molecular  inter¬ 
actions  . 
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SECTION  I 
INTRODUCTION 


The  viscosity  of  a  gas  in  the  low-density  limit  can  be 
described  by  the  Chapman -Enskog  theory  which  is  based  on 
Boltzmann's  equation  [1,2],  This  theory  is  akin  to  the 
perfect-gas  law  in  that  it  assumes  that  the  molecules  are 
randomly  distributed  in  configuration  space.  In  addition, 
the  theory  supposes  that  the  velocities  of  two  molecules 
which  are  about  to  collide  are  uncorrelated  (assumption  of 
molecular  chaos) . 

A  first  attempt  to  extend  the  theory  to  higher  densi¬ 
ties  was  made  by  D.  Enskog  for  the  case  of  a  gas  of  hard 
spheres  [3] .  In  this  theory  an  estimate  was  made  to  account 
for  the  effect  of  correlations  in  configuration  space?  these 
were  assumed  to  be  the  same  as  those  for  a  gas  in  equili¬ 
brium.  The  assumption  of  molecular  chaos  was  retained  for 
the  probability  distribution  in  velocity  space.  The  theory 
implies  the  existence  of  a  power  series  for  the  viscosity 
of  the  form 


n  =  no  +  nip  +  n2p2  +  n3p3  +  ...  (l) 

Such  a  series  is  analogous  to  the  virial  expansion  for  the 
equilibrium  properties  of  a  gas.  The  coefficients  of  the 
power  series  represent  the  effect  of  successively  higher- 
order  correlations  in  the  position  variables.  The  range  of 
these  correlations  is  closely  related  to  the  range  of  the 
molecular  interaction. 

A  proper  description  of  the  density-dependence  of  the 
transport  properties  also  requires  the  inclusion  of  an  esti¬ 
mate  of  the  deviations  from  molecular  chaos  in  the  velocity 
distribution.  Such  deviations  are  brought  about  by  sequen¬ 
ces  of  correlated  collisions  between  the  molecules  [4,5], 

In  the  past  it  was  naturally  assumed  that  these  effects 
could  also  be  represented  by  a  power  series  in  the  density. 
However,  the  range  of  such  velocity  correlations  is  not 
restricted  to  the  range  of  the  molecular  interaction,  but 
extends  to  distances  of  the  order  of  the  mean  free  path. 

Since  the  mean  free  path  is  itself  a  function  of  the  density, 
it  is  not  evident  that  a  power  series  which  is  formed  with 
density-independent  coefficients  does  exist.  In  fact,  a 
systematic  generalization  of  the  Boltzmann  equation  indi¬ 
cates  that  the  coefficient  of  the  quadratic  term  contains 
a  contribution  which  is  proportional  to  the  logarithm  of 
density,  so  that  we  should  write  [4,6] 

n  =  Ho  +  TUP  +  ri2  p2  £np  +  T.2P2  +  ...  (2) 
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as  a  replacement  for  equation  (1).  Nevertheless,  in  the 
past  experimental  results  on  transport  properties  were 
habitually  represented  by  a  power  series  like  that  in 
equation  (1) .  As  a  result  of  the  theoretical  developments 
alluded  to  before,  the  adequacy  of  a  power  series  cannot 
be  taken  for  granted  and  a  critical  re-examination  of  the 
experimentally  established  density  dependence  of  these 
properties  is  called  for. 

H.  J.  M.  Hanley,  R.  D.  McCarty  and  J.  V.  Sengers  [7] 
undertook  an  investigation  to  determine  whether  existing 
experimental  data  would  favor  the  theoretically  derived 
representation  (2)  in  preference  to  the  power  series  (1). 

A  detailed  account  of  this  work  was  presented  in  AEDC-TR- 
69-68  [4] .  Data  on  viscosity  which  were  available  at  the 
time  did  not  enable  them  to  perform  such  a  discrimination 
owing  to  inadequate  precision.  J.  T.  F.  Kao  and  R.  Kobayashi 
[8]  and  J.  W.  Gracki,  G.  P.  Flynn  and  J.  Ross  [9]  reached  a 
similar  conclusion.  On  the  other  hand,  the  more  extensive 
data  for  thermal  conductivity  [10,11]  seemed  to  favor  the 
inclusion  of  the  logarithmic  term.  This  conclusion  was  con¬ 
firmed  by  B.  Le  Neindre  and  coworkers  [11] . 

In  spite  of  this  agreement,  we  hesitated  to  accept  the 
conclusion  that  the  presence  of  the  logarithmic  term  was 
verified  experimentally.  First,  in  one  case,  the  precision 
of  the  measurements  was  only  one  of  ±1%.  Secondly,  the 
density  intervals  of  both  sets  of  measurements  were  quite 
large,  and  this  necessitated  the  inclusion  of  data  pertain¬ 
ing  to  excessively  high  densities  (800  Amagat  units)  in 
order  to  extract  the  required,  statistically  significant, 
information.  The  latter  point  is  particularly  important 
owing  to  the  circumstance  that  the  omission  of  higher-order 
terms  beyond  those  quoted  in  equations  (1)  and  (2)  may  not 
be  justified  in  such  an  extended  range. 

The  preceding  reasons  induced  us  to  undertake  new 
measurements  of  a  transport  property  and  to  arrange  them  so 
that  they  would  be  conducive  to  a  more  stringent  statistical 
analysis.  We  chose  to  measure  the  viscosity  of  nitrogen, 
argon  and  helium  in  an  oscillating-disk  instrument  [12] 
because  we  could  obtain  in  it  a  precision  of  a  few  parts  in 
ten  thousand.  Furthermore,  we  were  able  to  make  measure¬ 
ments  at  very  closely  spaced  intervals  of  pressure  and  thus 
to  provide  a  large  population  of  data  in  a  moderate  density 
range . 


2 


AEDC-TR-71  -190 


SECTION  II 

EXPERIMENTS 


Measurements  were  performed  on  nitrogen  (N2),  argon 
(Ar) ,  and  helium  (He)  at  a  nominal  temperature  of  25°C  and 
over  a  range  of  pressures,  as  follows: 


n2 

up 

to 

p  = 

105 

atm 

and 

p  =  0.12 

g/cm 

Ar 

up 

to 

p  = 

100 

atm 

and 

p  =  0.17 

g/cm 

He 

up 

to 

p  = 

100 

atm 

and 

p  =  0.02 

g/cm 

They  were  evaluated  on  a  relative  basis  [13,14].  The  tempe¬ 
rature  in  the  instrument  was  maintained  by  thermostating  the 
room  in  which  the  instrument  was  accommodated.  The  charac¬ 
teristics  of  the  suspension  system  are  given  in  Table  I. 

The  working  equation  [14] 


irpR^d 

21 


+  (26)'1/2 


(1  -  |a2> 


1  irpR2b3 
90  I 


Cxo  ' 


(3) 


1-3A2 
— TZ- 


0 


contains  the  calibration  factor  C  which  is  a  unique  function 
of  the  boundary- layer  thickness  [14] 


<5  =  VnT  /it  p  • 
o 


(3a) 


Here,  F  denotes  the  radius  of  the  disk,  d  is  its  thickness, 

I  is  its  moment  of  inertia,  bi  and  b2  are  the  upper  and 
lower  spacings,  respectively,  both  equal  (bi  =  b2  =  b) ,  T 
is  the  period  of  oscillation  in  vacua ,  0  =  T/T  is  the  0 
ratio  of  the  actual  period  to  that  in  a  vacuum,  and  A  is 
the  logarithmic  decrement  of  the  damped  harmonic  oscillation 
performed  by  the  system,  A0  denoting  its  value  in  vacuo. 
Finally 


x 


0 


(3b) 


3 


AEDC-TR-71-190 


TABLE  I 

CHARACTERISTICS  OF  SUSPENSION  SYSTEM 


Suspension  wire 


0.002  in.  diam. 
92%  Pt-8%  W. 
stress  relieved 


Wire  damping 

A  . 
o 

Total  separation  between 

n  = 

plates 

Upper  and  lower 

b  = 

separations 

Disk  radius 

R  = 

Disk  thickness 

d  = 

0.000040  ±  0.000004 

(0.28443  ±  0.00005)  cm 

bi  =  b2  =  (0.09006  ±  0.00005)  cm 

(3.4906  ±  0.001)  cm 

(0.10431  ±  0.00012)  cm 


Moment  of  inertia  of 
suspension  system 


I  =  (53.6032  ±  0.0006)  g*cm2 


Natural  periods  of 
oscillation 


Tq  =  (29.244  ±  0.002)  sec  at  25°C 
(first  calibration) 

To  =  (29.216  ±  0.002)  sec  at  25°C 
(second  calibration) 
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computed  numerically  from  equation  (3)  ,  determines  the 
viscosity,  n- 

In  order  to  perform  the  three  series  of  measurements, 
it  turned  out  to  be  necessary  to  calibrate  the  instrument 
twice;  owing  to  a  mishap,  the  suspension  wire  was  damaged 
before  the  measurements  on  argon  were  completed,  and  the 
oscillating  system  had  to  be  re-assembled.  The  points 
obtained  with  the  aid  of  the  second  calibration  are  marked 
with  asterisks  in  the  tables  of  results. 

The  calibration  curve  C(6)  was  determined  with  respect 
to  nitrogen  in  the  pressure  range  1-25  atm  and  at  25°C.  The 
calibration  values  were  taken  directly  from  Table  XV  of 
reference  [12] .  Since  the  temperature  in  the  instrument  was 
not  exactly  equal  to  25°C,  differing  from  it,  however,  by  at 
most  0.3°C,  a  small  correction  was  applied.  We  based  this 
correction  on  the  circumstance  that  the  excess  viscosity  is 
a  function  of  density, 

n(p,T)  -  no(0,T)  =  f(p)  (4a) 

in  both  representations  [151,  equations  (1)  and  (2).  Hence, 
along  an  isochore  p  =  const,  we  find  that 

n(p,T*)  =  n(p,T)  -  (rio(0,T)  -  r)  o  (0  ,T*)  } ,  (4b) 


where  T. =25°C.  In  this  manner  the  only  correction  that  needs 
to  be  applied  is  deduced  from  the  zero-density  values,  no, 
a  linear  estimate  of  0.47  yPoise/°C  being  adequate  for  the 
purpose.  The  corrected  value  now  refers  to  the  actual  density 
maintained  in  the  instrument  but  to  a  slightly  different 
temperature.  To  account  for  this  fact,  the  pressure  recorded 
during  the  measurement  must  be  corrected  accordingly. 

The  density  p,  was  calculated  from  the  equation 


=  1  +  BiP  +  B2P2  +  B3P3  (5) 

with 

R  =  82.082  atm  cm3/g  K 

M  =  28.014  g/mole  (on  carbon  12  scale). 
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The  virial  factors  Bi,  B2  ,  B3  listed  in  Table  II  were  inter¬ 
polated  from  NBS  Circular  564  [16] ,  and  a  further  linear 
interpolation  routine  between  20°C  and  30°C  was  programmed 
in  conjunction  with  the  working  equation  (3).  In  the  latter, 
the  values  of  C  were  evaluated  from  known  values  of  xQ,  the 
reverse  procedure  being  used  for  the  determination  of  the 
viscosity,  n, in  measurements  proper. 

The  two  calibration  curves  are  seen  plotted  in  Figure  1; 
they  can  be  correlated  as  follows: 


Ci  =  1.089  0645 


0.002  5248  0.000  7367  0.000  0334 

5  +  0.07  (6  +  0.07)  2  "  (6  +  0.07)  3 


(standard  deviation  1.5  x 


; 


C2  =  1.089  5202 


0.003  1036  0.000  5315  0.000  0217 

<$  +  0.0  7  (6  +  0.07)  2  "  (6  "+"  0.07)3  1 


(standard  deviation  1.8  * 


; 


The  experimental  results  are  given  in  Tables  III-V. 

The  density  of  nitrogen  was  evaluated  in  the  same  way  as 
during  a  calibration,  that  of  argon  and  helium  was  evaluated 
in  like  manner,  except  for  the  virial  coefficients  which  are 
also  listed  in  Table  II. 

The  data  for  argon  were  taken  from  the  measurements 
performed  by  E.  Whalley,  Y.  Lupien  and  W.  G.  Schneider  [17], 
those  for  helium  were  taken  from  the  work  of  W.  G.  Schneider  and 
J.  A.  H.  Duffie  [18] ;  it  proved  to  be  sufficient  to  retain 
terms  up  to  0 (P 2 )  only. 

In  all  cases  a  small  correction  for  temperature  was 
applied  to  reduce  the  data  to  25°C  as  mentioned  earlier  in 
connection  with  the  calibration  curves .  The  actual  tempera¬ 
ture  corrections  were  deduced  from  runs  at  slightly  different 
temperatures,  and  turned  out  to  have  the  following  values: 

Ar  0.71  uPoise/°C 

He  0.46  uPoise/°C 

In  this  manner,  the  standard  viscosity,  listed  in  the  last 
column  of  Tables  III-V,  represents  the  viscosity  of  the  gas  at 
25°C  and  at  the  actual  density  that  prevailed  during  the 
measurement.  For  this  reason,  we  evaluated  the  corresponding, 
corrected  pressure  and  listed  it  in  the  second  column  in  the 
tables,  as  already  mentioned. 
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TABLE  II  -  VIRIAL  COEFFICIENTS  FOR  NITROGEN ,  ARGON ,  AND  HELIUM 


Gas 

Temperature ,  °C 

Bi,  (atm)”1 

_2 

B2 ,  (atm) 

B3 ,  (atm)  3 

N  * 

20 

-2.42  x  10"4 

2.13  x  lo-6 

3.84  x  10”9 

30  ^ 

-1.47  x  10”4 

2.04  x  lo“6  1 

2.64  x  10”9 

Ar** 

20 

-7.5194  x  10”4 

2.389  x  10”6 

—  -  — 

30 

-6.2381  x  10”4 

2.1046  x  l(T6 

-  -  - 

He*** 

20 

4.633  x  10”4 

—  —  — 

—  —  — 

30 

4.471.  x  10"4 

-  -  - 

-  -  - 

From  reference  [16] 

** 

From  reference  [17] 

** 

From  reference  [18] 
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Graphs  of  standard  viscosity  versus  density  are  presented 
in  Figures  2,  3  and  4.  We  estimate  that  the  accuracy  of  the 
measurements  is  of  the  order  of  ±0.2%.  The  precision  and  in¬ 
ternal  consistency  is,  of  course,  much  better,  as  will  emerge 
from  the  analysis. 

In  assessing  the  reliability  of  the  preceding  experimental 
values,  it  is  necessary  to  realize,  as  pointed  out  in  our  ear¬ 
lier  experimental  work,  that  equation  (3)  is  quite  insensitive 
to  the  assumed  value  of  density.  This  is  due  to  the  fact  that 
the  damping  in  an  oscillating-disk  viscometer  which  is  equip¬ 
ped  with  closely  spaced  fixed  plates  is  very  nearly  independent 
of  density.  For  example,  if  the  density  in  equations  (3)  and 
(3a)  were  changed  by  as  much  as  1%,  the  viscosity  evaluated 
from  equation  (3b)  would  change  by  less  than  0.1%.  Thus,  an 
accurate  knowledge  of  the  density  is  inessential  for  the  eval¬ 
uation  of  viscosity.  By  contrast,  the  values  of  density  as¬ 
signed  to  the  corrected  pressures  at  25°C  must  be  computed 
with  great  care  when  an  analysis  of  the  density  expansion  of 
viscosity  is  undertaken.  Nevertheless,  the  analysis  is  not 
sensitive  to  uncertainties  in  the  actual  values  used  for  the 
virial  coefficients  in  equation  (5) .  A  choice  of  other 
values  for  these  virial  coeficients  would  represent  a  coordi¬ 
nate  transformation  p*  =  f(p)  where  f(p)  is  a  power  series  in 
p.  Such  a  coordinate  transformation  does  not  affect  the 
presence  or  absence  of  a  logarithmic  term  in  equation  (2) . 
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Density,  Amagat 


Figure  2.  Viscosity  of  nitrogen  as  a  function  of  density 


Poise 


Density,  Amagat 

Figure  4.  Viscosity  of  helium  as  a  function  of  density. 


TABLE  Ill 


THE  VISCOSITY  OF  NITROGEN 


Point 

number 

n 

Experi¬ 

mental 

pressure 

Pq , atm 

Corrected 

pressure 

P,atm 

Boundary 

layer 

thickness 

6 ,  cm 

Density 

p,g/cm3 

* 

Density 

p , Amagat 

Experi¬ 

mental 

temp. 

Te'°C 

Viscosity 

ne,yPoise 

Standard 

viscosity 

(at_25°C) 

n , pPoise 

1 

1.000 

1.000 

0.85037 

0.0011453 

0.91590 

24.920 

177. 94 

177. 98 

2 

1.714 

1.716 

0.64939 

0.0019645 

1.5710 

24.864 

178. 0Q 

178 .  0  c 

6 

3 

2.300 

2.302 

0.56057 

0.0026365 

2.1084 

24.727 

178. 0Q 

178. 13 

4 

3.009 

3.009 

0.49049 

0.0034470 

2.7566 

24.959 

178. 18 

178. 19 

5 

3.715 

3.718 

0.44124 

0.0042597 

3.4065 

24.722 

178. 19 

178 . 32 

6 

4.341 

4.344 

0.40438 

0.0049765 

3.9797 

24.819 

178. 32 

178. 4l 

7 

4.979 

4.981 

0.38142 

0.0057072 

4.5641 

24.890 

178. 3g 

178*44 

8 

5.579 

5.583 

0.36034 

0.0063979 

5.1164 

24.784 

178. 49 

178. 5g 

9 

6.613 

6.616 

0.33116 

0.0075827 

6.0639 

24.874 

178. 6? 

178. 73 

10 

7.307 

7.312 

0.31503 

0.0083814 

6.7027 

24.804 

178. 72 

178.8-^ 

11 

8.029 

8.036 

0.30053 

0.0092126 

7.3674 

24.744 

178. 7g 

178. 9q 

12 

8.954 

8.964 

0.28460 

0.010279 

8.2202 

24.651 

178. 8g 

179*°4 

13 

9.931 

9.922 

0.27082 

0.011378 

9.0991 

25.280 

179. 3q 

179. 1? 

14 

12.567 

12.579 

0.24056 

0.014433 

11.542 

24.702 

179. 45 

179. 5g 

15 

14.880: 

14.864 

0.22165 

0.017057 

13.641 

25.334 

180.0  . 

179. 8q 
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TABLE  III  -  THE  VISCOSITY  OF  NITROGEN  (continued) 


Point 

number 

n 

Experi¬ 

mental 

pressure 

Pe , atm 

Corrected 

pressure 

P,atm 

Boundary 

layer 

thickness 

6  fcm 

Density 

p,g/cm3 

* 

Density 

p ,Amagat 

Experi¬ 

mental 

temp. 

m  on 

e ' 

Viscosity 

ne,yPoise 

Standard 
viscosity 
(at  25°C) 

n,yPoise 

16 

20.051 

20.068 

0.19095 

0.023045 

18.429 

24.749 

180. 54 

180. 66 

17 

24.644 

24.638 

0.17274 

0.028307 

22.637 

25.068 

181. 48 

181. 45 

18 

29.838 

29.806 

0.15753 

0.034258 

27.396 

25.318 

182. 64 

182. 4g 

19 

34.816 

34.807 

0.14609 

0.040016 

32.001 

25.084 

183. 4g 

183. 45 

20 

39.919 

39.958 

0.13663 

0.045946 

36.743 

24.712 

184. 2g 

184. 43 

21 

44.682 

44.663 

0.12966 

0.051357 

41.070 

25.131 

185. 52 

185. 4g 

22 

49.853 

49.840 

0.12310 

0.057304 

45.826 

25.077 

186. 58 

186. 54 

23 

55.126 

55.106 

0.11750 

0.063344 

50.657 

25.113 

187. 8g 

187. 84 

24 

60.603 

60.602 

0.11241 

0 .069635 

55.688 

25.003 

189. 05 

189. 05 

25 

65.434 

65.419 

0.10862 

0.075132 

60.083 

25.070 

190. 44 

190. 4l 

26 

70.401 

70.421 

0.10504 

0.080825 

64.636 

24.916 

191. 6-l 

191. 65 

27 

73.599 

73.560 

0.10305 

0.084386 

67.484 

25.162 

192. 55 

192. 4? 

28 

80.403 

80.375 

0.09911 

0.092094 

73.648 

25.102 

194. 37 

194. 32 

29 

84.009 

83.970 

0.09727 

0.096140 

76.884 

25.139 

195. 44 

195. 3? 

30 

90.643 

90.624 

0.09414 

0.10360 

82.850 

25.062 

197.2, 

197. 2n 
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TABLE  III  -  THE  VISC0SIT3  OF  NITROGEN  (continued) 


Point 

number 

n 

Experi¬ 

mental 

pressure 

Pe,atm 

Corrected 

pressure 

P  ,atm 

Boundary 

layer 

thickness 

S  ,cm 

Density 
p, g/cm 3 

* 

Density 

p, Amagat 

Experi¬ 

mental 

temp. 

Te'°C 

Viscosity 

ne,uPoise 

Standard 
viscosity 
(at  25°C) 

n,pPoise 

31  ' 

94.555 

94.530 

0.09248 

0.10795 

86.328 

25.079 

198. 34 

198. 3q 

32 

97.821 

97.783 

0.09120 

0.11156 

89.215 

25.115 

199. 36 

199. 3L 

33 

105.237 

105.251 

0.08852 

0.11979 

95.797 

24.957 

201. 6C 

D 

201. 6? 

* 


1 


Amagat  =  0.00125046  g/cm3 
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TABLE  IV  -  THE  VISCOSITY  OF  ARGON 


Point 

number 

n 

Experi¬ 

mental 

pressure 

PQ,atm 

Corrected 

pressure 

P,atm 

Boundary 

layer 

thickness 

6  ,cm 

Density 

p,g/cm3 

Density** 

p ,Amagat 

Experi¬ 

mental 

temp. 

Te,‘C 

Viscosity 

nG,|iPoise 

Standard 
viscosity 
{at  25°C) 

n /pPoise 

1 

1.000 

1.000 

0.80283 

0.0016340 

0.91597 

24.911 

226. 28 

226. 34 

2 

*  1.497 

1.498 

0.65559 

0.0024472 

1.3718 

24.872 

226. 2q 

226. 2g 

3 

2.041 

2.041 

0.56212 

0.0033360 

1.8701 

25.040 

226. 48 

226. 45 

4 

*  2.480 

2.483 

0.50924 

0.0040593 

2.2755 

24.687 

226. 3g 

226 . 6^ 

5 

3.034 

3.034 

0.46106 

0.0049633 

2.7823 

24.982 

226. 6g 

226. 70 

6 

*  3.500 

3.496 

0.42949 

0.0057211 

3.2071 

25.302 

226. 95 

226. 74 

7 

4.048 

4.050 

0.39893 

0.0066301 

3.7166 

24.837 

226. 7Q 

226. 82 

8 

*  4.514 

4.509 

0.37817 

0.0073832 

4.1388 

25.315 

227. 0g 

226. 8g 

9 

5.110 

5.113 

0.35518 

0.0083756 

4.6951 

24.826 

227 . 0X 

227. 13 

10 

*  5.487 

5.484 

0.34287 

0.0089853 

5.0369 

25.155 

227. ly 

227. 0g 

11 

6.  035 

6.032 

0.32714 

0.0098873 

5.5425 

25.141 

227. 35 

227. 25 

12 

*  6.501 

6.496 

0.31508 

0.010651 

5.9706 

25.212 

227. 3g 

227. 24 

13 

7.158 

7.155 

0.011737 

6.5794 

25.124 

227. 5? 

227. 4g 

14 

*  7.488 

7.484 

0.29354 

0.012279 

6.8832 

25.150 

227. 53 

227. 42 

15 

*  8.498 

8.581 

0.27564 

7.8160 

25.177 

227 . 8- 

227. 6Q 

ost-it-moaav 


TABLE  IV 


Experi- 

Boundary 

Point 

mental 

Corrected 

layer 

number 

pressure 

pressure 

thickness 

n 

Pe,atm 

- - -  —  -  i 

P  ,atra 

6 ,  cm 

THE  VISCOSITY  OF  ARGON  (aontinued) 
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TABLE  IV  -  THE  VISCOSITY  OF  ARGOS  (continued) 


Point 

number 

Experi¬ 

mental 

pressure 

Corrected 

pressure 

Boundary 

layer 

thickness 

Density 

Density** 

Experi¬ 

mental 

temp. 

m  op 

e ' 

Viscosity 

ne,liPoise 

Standard 
viscosity 
(at  25°C) 

n#pPoise 

n 

P  ,atm 
e 

P,atm 

6  ,cm 

p.g/cm3 

p , Amagat 

31 

57.779 

57.725 

0.10741 

0.097371 

54.583 

25.280 

241.3. 

o 

241. 16 

32 

60.154 

60.122 

0.10535 

0.10152 

56.909 

25.161 

242. 06 

241. 95 

33 

65.162 

65.144 

0.10146 

0.11023 

61.792 

25.083 

243 .8 . 

4 

243. 7g 

34 

71.421 

71.380 

0.09729 

0.12108 

67.874 

25.172 

246. 23 

246. 1L 

35 

75.300 

75.297 

0.09491 

0.12791 

71.702 

25.012 

247. 56 

247. 55 

36 

79.926 

79.968 

0.09233 

0.13607 

76.277 

24.842 

249. 25 

249. 3g 

37 

85.029 

85.020 

0.08986 

0.14491 

81.232 

25.031 

251. 42 

251. 4q 

38 

90.303 

90.274 

0.08751 

0.15412 

86.395 

25.095 

253. 5? 

253. 5q 

39 

95.303 

95.253 

0.08548 

0.16285 

91.289 

25.156 

255. 6? 

255. 56 

40 

100.032 

10  0.084 

0.08365 

0.17132 

96.037 

24.845 

257. 55 

257. 66 

h  **  * 

Calculations  performed  using  calibration  2.  1  Amagat  =  0.0017839  g/cm3 
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TABLE  V  -  THE  VISCOSITY  OF  HELIUM 


Point 

number 

n 

Experi¬ 

mental 

pressure 

Pe , atm 

Corrected 

pressure 

P,atm 

Boundary 

layer 

thickness 

6  ,cm 

Density 

p,g/cm* 

* 

Density 

p ,Amagat 

Experi¬ 

mental 

temp. 

T  °c 

Viscosity 

T)e,yPoise 

Standard 
viscosity 
(at  25°C) 

n ,pPoise 

1 

1.190 

1.191 

2.1798 

0.00019465 

1.0905 

24.807 

198. 71 

198. 8q 

2 

2.014 

2.016 

1.6757 

0.00032937 

1.8452 

24.751 

l98-V 

198. 85 

3 

3.000 

3.003 

1.3729 

0.00049045 

2.7476 

24.718 

198. 

198. 74 

4 

4.021 

4.026 

1.1859 

0.00065719 

3.6817 

24.662 

198. 59 

198. 74 

5 

5.001 

5.007 

1.0639 

0.00081704 

4.5773 

24.647 

198 . 55 

198. 72 

6 

6.008 

6.015 

0.97070 

0.00098111 

5.4964 

24.644 

198. 63 

198. 7y 

7 

6.987 

6.993 

0.90039 

3.0011401 

6.3871 

24.734 

198. 5g 

198. 7X 

8 

8.028 

8.036 

0.83996 

0.0013096 

7.3367 

24.693 

198. 51 

198.6. 

O 

9 

9.049 

9.057 

0.79142 

3.0014752 

8.2644 

24.744 

198. 52 

198. 64 

10 

9.954 

9.965 

0.75460 

0.0016225 

9.0896 

24.658 

198. 5q 

198. 66 

11 

12.397 

12.404 

0.67682 

0.0020173 

11.301 

24.839 

198. 54 

198. 62 

12 

14.866 

14.867 

0.61852 

0.0024153 

13.531 

24.973 

198. 53 

198.5. 

4 

13 

17.227' 

17.230 

0.57487 

0.0027962 

15.665 

24.943 

198. 54 

198.5, 

0 

14 

19.622 

19.627 

0.53884 

0.0031816 

17.824 

24.930 

198. 48 

198. 51 

15 

24.882 

24.885 

0.47920 

0.0040244 

22.546 

24.966 

198. 5C 

198.5_ 
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TABLE  V  -  THE  VISCOSITY  OF  HELIUM  ( continued ) 


Point 

number 

n 

Experi¬ 

mental 

pressure 

P  ,  atm 
e 

Corrected 

pressure 

P,atm 

Boundary 

layer 

thickness 

6  ,cm 

Density 
p, g/cm 3 

* 

Density 

p , Amagat 

Experi¬ 

mental 

temp. 

Te,«C 

Viscosity 

r)e,pPoise 

Standard 
viscosity 
(at  25°C) 

H/WPoise 

16 

29. 965 

29.968 

0.43710 

0.0048353 

27.089 

24.975 

198. 4g 

!98.50 

17 

34.918 

34.921 

0.40543 

0.0056220 

31.496 

24.977 

198. 5C 

198. 56 

18 

39.919 

39. 918 

0.37948 

0.0064122 

35.923 

25.007 

198. 38 

198. 3g 

19 

45.090 

45.097 

0.35747 

0.0072274 

40.490 

24.955 

198. 44 

!98.46 

20 

50.002 

49.990 

0.33990 

0.0079940 

44.784 

25.073 

198. 44 

!98.40 

21 

55.160 

55.186 

0.32381 

0.0088047 

49.326 

24.862 

198. 35 

198. 4X 

22 

60.134 

60.148 

0.. 31062 

0.0095753 

53.643 

24.932 

198. 5Q 

198. 53 

23 

64.855 

64.889 

0.29920 

0.010308 

57.748 

24.845 

198. 2? 

198.3. 

4 

24 

69.857 

69.853 

0.28877 

0.011073 

62.034 

25.017 

198. 3? 

198. 3? 

25 

75.061 

75.109 

0.27876 

0.011878 

66.543 

24.811 

198. 32 

198. 4X 

26 

79.960 

79.963 

0.27052 

0.012619 

70.695 

24.987 

198. 4X 

198. 41 

27 

85.165 

85.207 

0.26233 

0.013416 

75.160 

24.854 

198. 35 

198. 42 

28 

90.064 

90.080 

0.25542 

0.014153 

79.289 

24.947 

198-38 

198*40 

29 

94.589 

94.712 

0.24934 

0.014851 

83.199 

24.613 

198. 3? 

198. 55 

30 

100.049 

100.152 

0.24271 

0.015667 

87.770 

24.692 

198. 2„ 

198.4, 

*  , 
1  Amagat  =  0.0001785  g/cm3 


AEDC-TR-71-190 


SECTION  III 


AEDC-TR-71-190 


METHOD  OF  ANALYSIS 


The  verification  of  a  proposed  mathematical  form  for 
the  density-dependence  of  the  viscosity  is  hampered  by  the 
fact  that  a  series  expansion  contains  an  infinite  number 
of  terms.  Moreover,  no  rigorous  expressions  are  available 
for  the  numerical  evaluation  of  the  coefficients  in  such 
an  expansion  beyond  the  Chapman-Enskog  zero-density  value 
of  no-  Thus  there  exists  no  a  priori  information  regarding 
the  number  of  terms  of  the  expansion  that  would  be  adequate 
in  representing  the  viscosity  in  a  specified  density  inter¬ 
val.  Even  the  dilute-gas  value,  n 0 ,  cannot  be  established 
independently  with  sufficient  precision,  and  its  most 
probable  value  must  also  be  extracted  from  the  experimental 
data. 


Nevertheless ,  it  is  possible  to  formulate  heuristic 
criteria  that  ought  to  be  satisfied,  if  a  postulated  equation 
is  to  be  consistent  with  a  population  of  experimental  data. 

In  Chapter  IV  of  AEDC-TR-69-68 ,  prepared  in  collaboration 
with  H.  J.  M.  Hanley  and  R.  D.  McCarty  [4,7],  such  a  set  of 
criteria  was  formulated  in  conjunction  with  a  systematic  al¬ 
gorithm  for  the  determination  of  successively  higher-order 
terms  in  an  equation  such  as  (1)  or  (2) ,  but  truncated  at 
some  point. 

In  order  to  understand  the  algorithm,  we  consider  the 
succession  of  polynomials 


no 

+ 

nip 

(7a) 

no 

+ 

nip 

+  n2P2 

(7b) 

no 

+ 

nip 

+  n2Pz  +  na  p3 

(7c) 

each  containing  one  more  term  than  its  predecessor.  We  now 
fit  the  first  equation  to  the  first  n=5  points  in  the  table 
of  results  for  a  particular  gas,  and  record  the  least-square 
values  of  no  and  m  together  with  their  standard  deviations 
and  the  variance,  a^,  for  the  set.  This  procedure  is  con¬ 
tinued  with  n  increased  by  unity  in  each  successive  fit.  It 
is  clear  that  at  first  the  variance,  an,  decreases  with' n, 
but  after  a  certain  value  n=n '  (or  p=p’)  has  been  reached, 
the  variance,  having  passed  through  a  stationary  value, 
begins  to  increase,  because  the  quadratic  term  comes  into 
play.  The  same  occurs  with  respect  to  the  standard  devia¬ 
tions  of  the  coefficients.  The  calculation  is  then  stopped, 
and  the  same  procedure  is  followed  with  equations  (7b)  and 
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(7c)  except  that  progressively  larger  values  of  n  and  n' 
are  used.  In  this  manner  the  algorithm  allows  us  to  deter¬ 
mine  an  interval  0  <  p  <  p'  which  provides  the  best  fit  of 
the  data  to  a  particular  equation.  As  a  last  step,  we 
repeat  the  same  calculation  in  relation  to  the  alternative 
form 


ti  =  no  +  niP  +  Ti2  p2  Anp  +  n2P2 


(7d) 


In  addition,  detailed  deviation  plots  are  prepared  to  ensure 
that  the  distribution  of  the  former  is  random,  as  expected. 

We  say  that  an  equation  is  consistent  with  the  experimental 
data  if  the  following,  obviously  necessary,  conditions  are 
satisfied : 

(i)  The  deviations  of  the  experimental  data  from 
those  calculated  with  a  fitted  equation  should 

be  random,  so  that  the  standard  deviation  can  be  interpreted 
as  true  measure  of  the  experimental  precision. 

(ii)  When  the  values  returned  for  the  coefficients  of 
the  equation  are  studied  as  a  function  of  n  (or  p) , 

they  should  be  independent  of  n  within  the  limits  prescribed 
by  their  respective  standard  deviations. 


(Hi)  Within  their  standard  deviation  the  coefficients 
should  not  change  when  the  next  term  is  added  to 
the  equation  and  the  density  interval  is  increased 
correspondingly . 


It  should  be  noted  that  the  algorithm  ensures  that 
criteria  (i)  and  (ii)  are  satisfied.  The  final  decision  as 
to  whether  a  power  series  or  the  series  with  the  logarithmic 
term  included  yields  a  superior  empirical  representation  is 
made  with  reference  to  criterion  (Hi)  . 
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SECTION  IV 


AEDC-TR-71  -190 


NUMERICAL  EXPERIMENTS 


Before  applying  the  procedure  described  in  the  preceding 
section  to  a  set  of  experimental  results,  it  was  thought  ad¬ 
visable  to  perform  a  preliminary  test  on  a  set  of  artificially 
generated  numbers  and  so  to  ascertain  whether  the  algorithm 
is  capable  of  producing  a  cle.ar  discriminant  for  the  presence 
of  the  logarithmic  term.  For  this  purpose  we  generated  data 
from  the  following  three  equations 

n  =  226.09  +  0 . 199p  +  1.5  x  10~3p2  -  1.8  x  10-6p3  '  (8a) 

n  =  226.09  +  0 . 199p  -  0.2  x  10“3p2Jlnp  +2  x  10-3p2  (8b) 

n  =  226.09  +  0 . 199p  -  2.0  x  10“3p2£np  +  2  x  10~3p2  (8c) 


The  coefficients  in  these  equations  have  been  chosen  to  cor¬ 
respond  closely  to  the  case  of  argon.  Numerical  data  were 
generated  at  the  densities  (in  Amagat  units)  listed  in  Table 
IV  for  argon  with  a  random  error  added  on  to  them;  the  latter 
was  taken  from  a  normal  distribution  with  a  variance  of 
a=0 . 07* . 

Equations  (7)  were  fitted  to  the  data  thus  generated  by  the 
method  of  least  squares.  We  used  the  Omnitab  computer  program  for 
statistical  analysis  which  uses  a  Gram-Schmidt  orthogonalization 
procedure  [19].  In  a  comparative  study  by  R.  H.  Wampler  [20]  this 
method  was  shown  to  be  reliable  for  the  purpose  at  hand. 

The  results  obtained  for  the  fitted  equations  were 
studied  as  a  function  of  p'  to  determine  the  range  0  <  p  <  p'; 
if  any,  in  which  criteria  (£)  and  (ii)  were  satisfied.  The 
most  probable  values  for  the  coefficients  of  equations  (7)  are 
presented  in  Table  VI.  In  this  and  the  subsequent  table  n 
refers  to  the  number  of  data  points. 

It  is  seen  that  the  data  generated  by  the  cubic  equation 
(8a)  produce  polynomials  (7a) ,  (7b)  and  (7c)  which  do  satisfy 
criterion  (Hi) .  More  precisely,  the  quadratic  equation 
reproduces  the  values  of  Tio  and  m  obtained  with  the  linear 
equation  and  the  cubic  equation  reproduces  the  values  of  no, 

Hi  and  n2  derived  from  the  quadratic  equation.  Thus  we 


The  authors  are  indebted  to  Dr.  D.  T.  Gillespie  for  providing 
them  with  suitable  sets  of  random  numbers. 
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TABLE  VI  -  TEST  OF  THE  ALGORITHM 


Data 

generated  n  p ' 

with 

note 

Ho 

T11±am 

(r)2±o_  )  .10s 

M  2 

(n3±an3) .10® 

(n i  ±o'  )  .io3 

02 

: 

°n 

12  6 

226.08+0.05 

0.202+0.014 

-  - 

-  - 

0.08 

30  52 

226.0710.03 

0.20210.003 

1.3710.07 

-  - 

-  - 

0.07 

Eqn.  (8a)  40  96 

226.0910.03 

0.19810.004 

1.510.1 

-2.110.7 

-  - 

0.07 

40  96 

226.10+0.03 

0.19210.006 

2.510.4 

-  - 

-0.2310. 09 

0.08 

(exact) 

(226.09) 

(0.199) 

(1.5) 

(-1-8) 

(0) 

(0.07) 

13  7 

226.0810.05 

0.204+0.012 

_  _ 

—  _ 

0.07 

35  72 

226.0610.02 

0.20710.002 

1.0610.03 

-  - 

-  - 

0.07 

Eqn.  (8b)  40  96 

226.0710.03 

0.20410.004 

1. 2+0.1 

-2.010.7 

-  - 

0.07 

40  96 

226.09+0.03 

0.19810.006 

2.110.4 

-  - 

-0.2310.08 

0.07 

(exact) 

(226.09) 

(0.199) 

(2) 

(0) 

(-0.2) 

(0.07) 

15  8 

226.10+0.04 

0.18710.009 

-  - 

-  - 

-  - 

0.07 

24  28 

226.00+0.04 

0.24210.008 

-6.110.3 

-  - 

-  - 

0.08 

Eqn.  (8c)  31  55 

226.0310.04 

0.23010.008 

-4.810.4 

-3315 

-  - 

0.07 

40  96 

226.09+0.03 

0.19810.006 

2.110.4 

-  - 

-2.0310.08 

0.07 

(exact) 

(226.09) 

(0.199) 

(2) 

(0) 

(-2) 

(0.07) 

AEDC-TR-71-190 


AEDC-TR-7 1-190 


conclude  that  the  cubic  polynomial  is  consistent  with  the 
data  set.  The  value  m  =  0.192±0.006,  obtained  from  the 
logarithmic  equation  {7c)  ,  is  somewhat  lower  than  the  cor¬ 
responding  term  m=0.202±0.003  in  the  linear  equation. 
Nevertheless,  the  two  values  agree  within  their  combined 
standard  deviations  so  that  our  criteria  do  not  permit  us 
to  rule  out  a  logarithmic  term  when  it  is  multiplied  with 
a  coefficient  of  2  x  10“^.  This  conclusion  is  confirmed 
by  the  results  obtained  from  the  data  generated  according 
to  equation  (8b) .  This  equation  contains  a  logarithmic 
term  with  a  coefficient  of  this  magnitude.  According  to 
our  criteria,  both  the  cubic  equation  (7c)  and  the  logarith¬ 
mic  equation  (7d)  are  consistent  with  this  data  set. 

Finally,  we  considered  the  data  generated  by  equation 
(8c)  where  the  coefficient  of  p2Jinp  is  ten  times  larger 
than  that  in  the  previous  equation.  Our  procedure  now 
shows  very  clearly  that  the  polynomials  are  not  consistent 
with  the  data  set.  In  other  words,  the  quadratic  and  cubic 
equations  do  not  reproduce  the  values  of  no  and  m  of  the 
linear  equation,  nor  does  the  cubic  equation  reproduce  the 
value  of  n 2  of  the  quadratic  equation.  On  the  other  hand, 
the  logarithmic  equation  (7d)  does  reproduce  the  factors  no 
and  m  within  their  combined  standard  deviations  and  is  thus 
consistent  with  the  data. 

We  conclude  that  our  algorithm  can,  indeed,  distinguish 
between  a  polynomial  such  as  equation  (7c)  or  a  logarithmic 
equation  (7d) .  Of  course,  for  data  with  a  finite  precision, 
the  distinction  can  only  be  made  with  a  finite  resolution. 
With  our  current  relative  precision  of  3  x  10”^  we  can  no 
longer  detect  a  logarithmic  contribution  if  its  coefficient 
is  as  small  as  2  x  10”4. 
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SECTION  V 


ANALYSIS  OF  THE  EXPERIMENTAL  DATA 


The  experimental  viscosity  data  were  subjected  to  the 
analysis  described  in  the  previous  two  sections.  The  most 
probable  values  thus  deduced  for  the  coefficients  of  the 
equation  (7a),  (7b),  (7c)  and  (7d)  are  listed  in  Table  VII. 

In  fitting  the  equations  equal  weights  were  assigned  to  the 
individual  data  points.  However,  essentially  the  same 
results  are  obtained  when  a  weight  proportional  to  the 
value  of  viscosity  is  assigned  to  the  data,  as  shown  in 
Table  VIII.  We  note  that  the  analysis  indicates  that  the 
relative  precision  of  the  experimental  data  is  better  than 
3  x  10~4.  Deviation  plots  showing  nexp  ”  ^calc  for  nitrogen 
and  argon  are  presented  in  Figures  5,  6,  7,  8.  In  particu¬ 
lar  from  Figures  6  and  8  we  note  that  the  deviation  plots 
obtained  with  the  cubic  equation  (7c)  and  the  logarithmic 
equation  (7d)  are  virtually  indistinguishable  for  nitrogen 
as  well  as  for  argon.  Evidently  the  two  equations  satisfy 
criterion  ( i )  equally  well. 

If  we  apply  in  addition  criterion  (Hi) ,  we  see  from 
Table  VII  and  VIII  that  the  power  series  is  consistent  with 
the  experimental  data  for  nitrogen  and  argon  to  a  high  degree 
of  satisfaction.  This  means  that  the  linear,  quadratic  and 
cubic  equations  lead  to  the  same  coefficients.  The  logarith¬ 
mic  equation  (7d)  does  not  reproduce  the  value  of  m  of  the 
linear  equation  nearly  as  well.  Thus  we  conclude  that  from 
viscosity  data  with  a  relative  precision  of  3  x  lo-4  a 
p2£np  term  cannot  be  detected  in  the  case  of  nitrogen  or 
argon.  It  turns  out  that  in  the  case  of  helium  we  cannot 
obtain  significant  information  for  all  coefficients,  but  the 
results  for  the  linear  and  quadratic  equation  are  again  in 
satisfactory  agreement. 

As  demonstrated  in  the  previous  section  our  procedure 
would  have  revealed  the  existence  of  a  logarithmic  term  if 
its  coefficient  were  significantly  larger  than  2  x  10"4 
UPoise/Am2.  Therefore,  we  conclude  that  such  a  coefficient 
is  not  larger  than  the  above  value,  if  it  exists  at  all. 

In  order  to  make  this  statement  more  precise,  we  rewrite 
equation  (7d)  as 

f(n)  =  n  -  n2P2£np  =  no  +  tup  +  n 2 p 2  (9) 

In  this  new  equation  ri2  is  not  considered  as  a  coefficient 
to  be  determined,  but  as  a  parameter.  Each  value  for  the 
parameter  n 2  defines  a  new  set  of  data  points  f(n).  For  the 
"correct"  value  of  n 2  the  data  f(n)  should  satisfy  our 
criteria  when  fitted  to  a  polynomial.  Even  when  we  relax 
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our  criteria  to  allow  variations  within  twice  the  standard 
deviations,  we  conclude  that  values  for  n 2  such  that 
| n  2 |  >  6  x  10”4yPoise/Am2  cannot  possibly  be  reconciled  with 
our  experimental  data  for  nitrogen  or  argon. 
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TABLE  VII 


Statistical  analysis  of  the  experimental  data  (equal  weights) . 


n 

P ' 

no  * 

no 

111  *  V 

(02  1  0  ) .103 

nz 

(n,  ±  °n3>-106 

(02  *  on O .103 

°0 

Am 

yPoise 

yPoise/Am 

yPoise/Am2 

yPoise/Am3 

yPoise/Am2 

yPoise 

r 7 

5 

177 . 85*  0 . 02, 

0 . 133*  0 . 007 

-  - 

-  - 

-  - 

0.02 

21 

41 

177 . 85*  0 . 02 

0 . 1321 0.002 

1.31*0.06 

-  - 

-  - 

0.03 

33 

96 

177 . 85*  0 . 02 

0 . 131±  0 . 0 02 

1.36+0.06 

-1.5+0. 4 

-  - 

0.04 

'33 

96 

177.87*0.02 

0.126*0.004 

2. 1+0. 3 

-  - 

-0.18*0.04 

0.04 

■ 

f14 

7 

226 . 10*  0.04 

0 . 199*  0 .009 

M  __ 

_  _ 

—  _ 

0.06 

Ar< 

l27 

42 

226 .09*  0.02 

0 . 196*  0 .003 

1.59*0.09 

-  - 

-  - 

0.05 

140 

96 

226 . 08*  0.02 

0 . 198*  0.003 

1.54+0.09 

-1 . 9±  0 . 6 

-  - 

0.06 

■ 

^40 

96 

226. 09* 0.03 

0 . 192*  0.005 

2. 4+0. 4 

-  - 

-0 . 22±  0 .07 

0.06 

He< 

f14 

18 

198.83*0.02 

-0 . 0191 0 .002 

-  - 

-  - 

-  - 

0.03 

[22 

54 

198.82-0.02 

-0 . 0191 0.003 

0.24*0.05 

-  - 

* 

0.05 
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TABLE  VIII 


STATISTICAL  ANALYSIS  OF  THE  EXPERIMENTAL  DATA  (proportional  weights) 


no  ±  cr 


ni  ±  a. 


±  a  ) .103 

T1 2 

I  C  Q  /  A  TT»  ^ 


Am 

yPoise 

yPoise/Am 

yPoise/Am 

1  7 

5 

177.85+0.02 

0.133+0.007 

-  - 

21 

41 

177.85+0.02 

0.132+0.002 

1.31+0.06 

33 

96 

177 . 85±0 .02 

0.131+0.002 

1.37+0.06 

■33 

96 

177.87+0.02 

0.126+0.004 

2. 1+0.3 

14 

7 

226.10+0.04 

0.199+0.009 

. 

27 

42 

226.09+0.02 

0.196+0.003 

1.59+0.09 

40 

96 

226.08+0.02 

0 . 199±0 . 003 

1.54+0.09 

40 

96 

226.09+0.03 

0.192+0.005 

2. 4  +  0. 4 

14 

18 

198.83+0.02 

-0.019±0 .002 

-  - 

22 

54 

198 . 82±0 .02 

-0 . 019±0 .003 

0 . 24±0 . 05 

(ns  ±  an  ) .106 

'1  3 

yPoise/Am3 


-1 . 6±0 .4 


-1.9±0 .6 


(n2  ±  a  i) .103 

1)2 

yPoise/Am2 


-0.18±0.05 


-0.22+0.07 


0.01 

0.02 

0.02 

0.02 

0.03 

0.02 

0.03 

0.03 

0.02 

0.03 
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Figure  6.  Deviation  plots  for  the  viscosity  of  nitrogen 
(a)  Equation  (7c) ;  (b)  Equation  (7d) . 
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Fi.gu.fiz  7.  Deviation  plots  for  the  viscosity  of  argon 
(a)  Equation  (7a);  (b)  Equation  (7b). 
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Figure  8.  Deviation  plots  for  the  viscosity  of  argon 
(a)  Equation  (7c) ;  (b)  Equation  (7d) . 
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SECTION  VI 
DISCUSSION 


We  have  demonstrated  that  the  experimentally  determined 
density  dependence  of  the  viscosity  of  argon  and  nitrogen 
at  25°C  can  be  consistently  represented  by  a  power  series  in 
the  density.  Implicit  in  this  conclusion  is  the  assumption 
that  the  first  four  terms  of  the  density  expansion  (1)  and 
(2)  are  adequate  to  describe  the  data  in  the  density  range 
under  consideration.  The  existence  of  a  logarithmic  term 
of  the  form  p2£np  cannot  be  detected  on  the  basis  of  vis¬ 
cosity  data  whose  relative  precision  is  3  x  10“  .  If  such 
a  term  exists,  it  is  concluded  that  the  absolute  value  of 
its  factor,  ri2  ,  must  be  appreciably  smaller  than  6  *  10“4 
yPoise/Am2  for  argon  as  well  as  nitrogen  at  25°C.  From  a 
practical  point  of  view,  there  is  no  need  to  include  such  a 
term  in  describing  the  density-dependence  of  the  viscosity 
within  the  best  precision  currently  available . 

We  would  like  to  point  out  that  the  same  conclusion  is 
suggested  by  the  theoretical  results  recently  reported  in 
AEDC-TR-71-51 [5] .  Although  the  theory  does  predict  the 
existence  of  a  logarithmic  term  p2fi,np,  for  a  comparison  with 
experiment  it  is  important  to  have  an  estimate  of  its  factor, 
n2,  relative  to  the  magnitude  of  the  coefficient  ri2  of  the 
quadratic  term.  A  theoretically  calculated  value  for  ri2  is 
not  available  at  the  present.  Nevertheless,  we  venture  to 
make  an  estimate  based  on  some  considerations  developed  in 
AEDC-TR-71-51.  In  that  report  a  distinction  was  made  between 
correlations  that  are  of  a  statistical  nature  versus  those 
that  are  of  a  purely  dynamical  nature.  Statistical  cor¬ 
relations  are  related  to  those  molecular  configurations 
where  the  particles  are  close  together  in  their  mutual  inter¬ 
action  range.  The  dynamical  correlations  are  brought  about 
by  sequences  of  successive  binary  collisions  that  are 
spatially  separated  by  distances  larger  than  the  molecular 
interaction  range.  In  AEDC-TR-71-51  we  indicated  that  the 
contributions  from  purely  dynamical  correlations  are  con¬ 
siderably  smaller  than  those  from  statistical  correlations. 

It  was  shown  that  for  the  first  density  coefficient  m  the 
dynamical  correlations  constitute  at  most  a  few  percent 
relative  to  the  statistical  correlations,  referred  to  as 
double-overlap  term  or  Enskog  term. 

It  should  be  remembered  that  the  origin  of  the  p2&np 
term  in  (2)  is  entirely  of  a  dynamical  nature,  whereas  the 
coefficient  n 2  of  the  quadratic  term  will  be  determined 
primarily  by  statistical  correlations .  Thus  we  expect  that 
the  coefficient  n 2  would  be  considerably  smaller  than  the 
corresponding  coefficient  n  2  -  This  explanation  is  supported 
by  some  calculations  for  a  two-dimensional  gas  reported  in 
AEDC-TR-69-68  [4,21].  The  corresponding  ratio  for  a 
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two-dimensional  gas  m/ni  turned  out  to  be  as  small  as 
3  x  10“2,  we  think  that  this  result  is  indicative,  not¬ 
withstanding  the  fact  that  recent  developments  have  indi¬ 
cated  that  transport  coefficients  of  a  two-dimensional  gas 
may  not  exist  in  the 'long  wave  length 'limit  [22] . 

B^.sed  on  the  foregoing  considerations  we  estimate 
that  n 2 /n a  should  be  not  larger  than  a  few  percent.  Since 
ri2:slO-3yPoise/Am2  for  argon  and  nitrogen,  we  estimate  that 
H2  would  be  at  most  of  the  order  of  10“5yPoise/Amz .  This 
estimate  is  not  only  consistent  with  our  experimental  results, 
but  we  conclude,  in  retrospect,  that  a  logarithmic  term 
would  not  reveal  itself  even  if  the  precision  of  the  experi¬ 
ment  would  be  increased  by  one  or  two  orders  of  magnitude. 

Both  experiment  and  theory  indicate  that  a  power  series 
representative  for  the  density  dependence  of  the  transport 
properties  will  be  adequate  for  practical  purposes. 
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APPENDIX 


In  addition  to  the  principal  conclusions  of  the 
present  paper  regarding  the  form  of  the  density  expansion, 
the  experimental  data  can  be  utilized  in  a  practical  way 
to  provide  us  with  very  precise  data  on  the  zero-density 
value  of  the  viscosity  of  the  three  gases  and  on  the 
linear  correction  term.  Referring  to  Table  VII  we  can 
summarize  them  as  follows: 


NITROGEN: 

T)0 

= 

(177.85+0.02) 

yPoise  at  25°C 

Til 

= 

(0.132+0.007) 

yPoise/Amagat  at 

25°C 

ARGON: 

no 

= 

(226 . 09±0 . 04) 

yPoise  at  25°C 

m 

= 

(0.198+0.009) 

yPoise/Amagat  at 

25°C 

HELIUM: 

no 

= 

(198.82s±0.02) 

yPoise  at  25°C 

m 

= 

-(0.019±0.003) 

yPoise/Amagat  at 

25°( 

Evidently,  the  preceding  optimum  values  disregard  the  pos¬ 
sibility  of  the  existence  of  unaccounted  for  systematic 
errors . 
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